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We study the local magnetization in the 2-D Ising model at its critical temperature on a semi- 
infinite cylinder geometry, and with a nonzero magnetic field h applied at the circular boundary 
of circumference f3. This model is equivalent to the semi-infinite quantum critical 1-D transverse 
C^l ■ field Ising model at temperature T oc /3 _1 , with a symmetry-breaking field oc h applied at the point 

boundary. Using conformal field theory methods we obtain the full scaling function for the local 
magnetization analytically in the continuum limit, thereby refining the previous results of Leclair, 
<n: Lesage and Saleur in Ref. 1. The validity of our result as the continuum limit of the 1-D lattice 

■ model is confirmed numerically, exploiting a modified Jordan- Wigner representation. Applications 

' of the result are discussed. 

m 

7— I ■ 

I. INTRODUCTION 

^ ■ The Ising model is a classic paradigm of statistical mechanics, and continues to find powerful application in diverse 
areas of modern physics. 2 ~ 4 It also reveals unique and generic universal behavior associated with boundaries. 5-7 In its 
quantum 1-D chain version, the critical boundary Ising model (BIM) reads 
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H = - X>X+i + of] - h B <j z . (1) 
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The uniform field along x is fixed such that the bulk system is at the critical point between Ising order and the 
disordered phase. The symmetry a z <-> — a 7 " is broken when a finite magnetic field hs 7^ is applied at the point 
boundary. Such a boundary field cannot lead to a finite bulk magnetization. Importantly however, it does cause a 
renormalization group flow from a free boundary condition hs = to a fixed boundary condition Hb — > ±00. 

The renormalization group flow associated with this BIM has been shown to be at the heart of boundary critical 
phenomena occurring in a surprising variety of low-dimensional correlated electron systems, such as Luttinger liquids 
containing an impurity, 1 coupled bulk and edge states in non-abelian fractional quantum Hall states, 8,9 as well as 
quantum dots near the two impurity Kondo 10,11 or the two channel Kondo 12,13 critical points. 

In the continuum, the BIM is in fact integrable, 14 both in the massless bulk critical case, and also in the massive 
regime away from the critical point. Certain correlation functions can then be calculated exactly using Form Factor 
. methods; 15,16 although in the bulk critical case relevant to Eq. (1), many important quantities cannot be easily ob- 
■ tained due to the proliferation of many-particle excitations. On the other hand, Chatterjee and Zamolodchikov 17 (CZ) 
showed that conformal field theory imposes linear differential equations which fully determine correlation functions in 
this limit. Of course, conformal field theory has been used previously for systems with conformal-invariant boundary 
I conditions. 6,7 The remarkable feature of the result of CZ is that the correlation functions are still determined by 
7-H ■ differential equations even for non-conformal invariant boundary conditions obtained at finite boundary field. 

The method of CZ was applied to the calculation of magnetization as a function of distance from the boundary, x. 
. ,-h . On the semi-infinite plane, equivalent to the quantum 1-D model, Eq. (1), at zero temperature, their result reads 17 

Mx)) h = 2 13 / 8 ^hx 3/8 ^{l/2, 1; 8tt/i 2 x), (2) 

where ^ is a degenerate hypergeometric function. Here (o~(x)} has the standard field-theory normalization, which we 
emphasize is only proportional to (er|) of a particular lattice model, such as Eq. (1). Indeed, x oc j, and h oc hs 
provided that hs -C 1. At short distance one thus obtains, 

(a(x)) = -2 13 / 8 /ix 3/8 [ln(x) + 0(1)] + 0(x 7/s ) for x < 1, (3) 

and at long distances (cr(x)) (2/rr) 1 / 8 , corresponding exactly to the result for fixed boundary condition, obtained 
from boundary conformal field theory. 18 As such, the exact function, Eq. (2), captures the full crossover behavior 
between two fixed points where conformal invariant boundary conditions do hold. 

Since the problem for finite h does not in general possess conformal invariance at the boundary, it is not possible to 
generalize Eq. (2) to other geometries by means of a simple conformal mapping. However the method of CZ can be 
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applied directly to other geometries, yielding a new set of differential equations (this was demonstrated explicitly for 
the 2-D disk geometry by CZ 17 ). Similarly, Leclair, Lesage and Saleur 1 (LLS) applied the method to the geometry of 
a semi-infinite cylinder. In the present paper we shall be concerned with this semi-infinite cylinder geometry, whose 
boundary consists of a circle with circumference (3, at which the boundary field h is applied. This classical 2-D Ising 
model is equivalent to the quantum chain model Eq. (1) at finite temperature T oc /3 _1 (see e.g. Ref.19). 

We re-examine the result of LLS for the local magnetization in Sec. II. Whereas those nonperturbative results give 
the full x-dependence of the local magnetization for any h and /3, we find that a more general ansatz for the local 
magnetization allows for an additional multiplicative factor f(2/3h 2 ). The physical meaning of this missing factor is 
then explained. The full scaling function for the local magnetization is determined in Sec. Ill; while the lattice model 
Eq. (1) is studied directly in Sec. IV. The local magnetization on the lattice is calculated numerically, and the results 
compared with the refined analytic solution, showing excellent agreement. The paper ends with a short summary, 
where implications and applications of the results are discussed. 



II. REFINEMENT OF EARLIER RESULTS 



LLS considered a classical Ising model on the half-cylinder in the continuum limit. They calculated the local 
magnetization as a function of the distance x from the circular boundary of circumference /3, which was conveniently 
written in the form 1 




= -T-F7- 9(X), (4) 



with X = (1 — coth=|p)/2 and where (cr(x)) is independent of r e (0,/3) due to translation symmetry along the 
boundary. LLS derived a linear differential equation for g(X), which reads 1 



(dX) 

parametrized in terms of A = 2/3 h 2 . Their solution is 1 
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( X - X )7J^2 + [ 1 + -- 2X )—~7)9(X)=0, (5) 
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where 2-Fi(a, b, c, z) is the Gauss hypergeometric function. Below we will use its integral representation 

t\c\ r 1 t b - 1 (i-t) c - b - 1 

where T(y) is the gamma function. This result is normalized with an overall constant such that at long distances 
(<t(x)) — > ( sinh |9 i !J x ) recovers the expected result for fixed boundary conditions (taking /3 — > 00 then yields (cr(x)) — > 

(2/x) 1/s , consistent with Ref. 18). 

In this paper we point out that the differential equation Eq. (5) leaves a freedom which goes beyond an overall 
normalization constant. Unlike the zero-temperature case (corresponding to the semi-infinite plane, /3 —¥ 00), here 
the normalization of Eq. (6) can itself be a scaling function of A. We thus replace Eq. (6) by the more general ansatz, 

11 1 — coth \ 
( -,-;l + 2/3/ 1 2 , 2 _ ^~)' (8) 

which depends explicitly on the function /(2/3/i 2 ), determined in Sec. Ill, below. Eq. (8) implies /(oo) = 1, so that 
{a(x))h,fs recovers asymptotically the behavior of the fixed boundary condition fixed point. 

We note that the same subtlety occurs with other geometries, as highlighted by CZ in the case of the disk. 17 In that 
case, the additional scale in the problem is the disk radius R; and an additional scaling function of Rh 2 (analogous 
to our A = 2/3ft, 2 ) appears in the expression for the local magnetization. As in the present case, this function is not 
fixed by the linear differential equations. 17 
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Finally, we comment briefly upon the physical significance of the scaling function f(2(3h 2 ). It describes the depen- 
dence on the additional thermal scale influencing the renormalization group flow at T ^ 0. In accord with physical 
expectation, the renormalization group flow is cut off at the external scale given by max{T, x -1 }. Since h grows under 
renormalization (and has scaling dimension 1/2), 6,7 we should consider two regimes depending on the ratio between 
this external scale and the field-induced scale ~ h 



2. 



max{T, x 1 } h 2 : free boundary condition, 

max{T, a; 1 } <C h 2 : fixed boundary condition. (9) 

These regimes are illustrated in Fig. 1. The important consequence following from this is that at finite temperatures, 
the fixed boundary condition fixed point is not always reached on taking x — > oo. The single scaling function 
f(2(5h 2 ) thus describes the crossover from free to fixed boundary condition at a given x, upon decreasing temperature. 
Obviously its effect is most apparent at large x, since there is no crossover at small x. However, as suggested by 
Fig. 1, the system is always close to the free boundary condition fixed point at small x, and this fact will prove useful 
in determining f(2f3h 2 ), as considered in the next section. 
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Figure 1: Schematic phase diagram of the BIM as a function of temperature and distance from the boundary. The dashed line 
denotes the crossover between free and fixed boundary conditions occurring when max{T, x -1 } ~ h 2 . 



III. DETERMINATION OF THE SCALING FUNCTION /(2/3/i 2 ) 

In this section we find the function f(2/3h 2 ) appearing in Eq. (8). Since this function is a scaling function of j3h 2 
and does not depend on distance x, it could in principle be determined at any given x. While its influence is most 
pronounced at large x, where the system undergoes a crossover as function of T (see Fig. 1), here we determine 
f(2/3h 2 ) by exploring the small x behavior, where the system remains close to free boundary condition fixed point. 
Importantly, the resulting behavior at small x is perturbative in h regardless of fiti 2 , as shown explicitly below. 

First we note that at both large and small A, the short-distance behavior of (<r(x)) is linear in h. As A — > oo, one 
sees this directly from the small x expansion of the exact T = result of CZ, Eq. (3). In the opposite limit A —J- 0, the 
behavior is by definition perturbative in h, and so the leading correction to magnetization is of course also linear in 
h. In the next subsection, we perform first-order perturbation theory in the boundary field h, with respect to the free 
boundary condition fixed point. The key point is that its short-distance behavior yields precisely Eq. (3), implying 
that 

(<j{x)) h ^^-2 13 ^hx 3 ^[\n(x)+0{l)]+0(x 7 / s ) for x « /3, h~ 2 (10) 

holds at short distances x -C j3, hr 2 for any A. Naively one might think that the coefficient of the x 3 ^ 8 [ln(x) + 0(1)] 
term could be renormalized by higher orders in h. But the scaling form of the problem implies that every power of h 
is accompanied by a power of \ffi [or y/x which gives a subleading x dependence to Eq. (10)]. Such terms diverge as 
T — > 0, and so this renormalization is not consistent with the exact nonperturbative T = result, Eq. (2), which is 
well-behaved at short distances, Eq. (3). 

Finally, we consider the short distance expansion of our ansatz Eq. (8), which using Eq. (7) is found to be 

(a(x)) h , = -f(2/3h 2 )^ + g| x W Mx) + 0{1)] + 0{x 7/ 8) for x ^ i (n) 
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Comparing Eqs. (10) and (11) we now obtain the scaling function 



This function increases monotonically as shown in Fig 2, and has asymptotic limits /(A <C 1) 



(12) 

/ 7rA and /(A > 
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0.5 1. 1.5 2. 2.5 3. 
Figure 2: Plot of Eq. (12). 



We note that a similar perturbative method was used by CZ to fix the scaling function of Rh 2 for the disk geometry. 
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A. Perturbation theory in the boundary field h 

In this subsection we show that the form of Eq. (10) indeed follows from perturbation theory around the free 
boundary condition fixed point. We obtain the full x//3 dependence of the magnetization at small h, recovering 
perturbatively the A^O limit of Eq. (8). 

The continuum limit of the critical classical 2-D Ising model is described by a c = 1/2 conformal field theory, which 
admits a Lagrangian formulation in terms of the free massless Majorana Fermi field (ip,ip), with the action 



Sa = h I d2z[% 



(13) 



Here {z, z) = (t + ix, t — ix) are complex coordinates and d 2 z = drdx. In the presence of a boundary B with a 
magnetic field h, the action can be decomposed into a bulk part and a boundary part, 



S = 



1 

2tt~ 



d 2 z[^dgi; + i;d z ip} + h 



OB- 



(14) 



The boundary operator ob was identified in Refs. 6,7 with a dimension 1/2 operator ctb(t) ~ "0(7", x = 0), associated 
with the fermion field at the boundary. In our case B = dT> is a circle parametrized by t € [0,/?] and V is the 
semi-infinite cylinder. Following Cardy's method of images 6 the one point function of the magnetization is a(z±, z<i) = 
(a l{zi)o l{z2)) , where ol{z) is a dimension 1/16 left moving field living in the geometry of the infinite cylinder. We 
then obtain conformal invariant boundary conditions, with the 'boundary' at x = 0. The boundary field h is now 
considered as a perturbation to the free boundary condition fixed point. To first order in h, 



a^( Zl ,z 2 ) = h / dT(o( Zl )oB(0,T)o(z 2 )). 



(15) 



The 3-point function appearing in the integrand is fully determined by conformal invariance, and one obtains up to 
a normalization constant N 



cr (1) («i,«a) = hN 



sin j(z\ - z 2 ) 



3/8 



dr 



sin(f(r-zi))sin^|(r-^ 2 )) 



1/2' 



(16) 
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The physical magnetization is obtained by setting (cr(x)) — (<r(zi = ix, z-i = —ix)). We now take z\ = ix, z-i = —ix 
in Eq. 16 and use the trigonometric identity, 

2 sin (j(r - *)) sin g(r - z 2 )J = cos - cos j . (17) 



The integral in Eq. (16) then becomes 



L 



2lx df) 2?r / 1 1 , 1 - coth ' 



jl/2 +'"' 1/2 -cosg Jsinh^ V 2 2 



in terms of w = e 4lra: /^ and 9 = 27tt//3. The constant N was carefully accounted for by CZ. 17 Using this and 
Eq. (18), first-order perturbation theory in the boundary field h yields 

_/ \ 1/8 A i i_ co th^\ 

("W) ? = h ^ l^f) 2Fl [r 2 ; 1; 2-^ J + ° {h2) - (19) 

It is interesting to compare this with the full result of LLS, Eq. (6). At small h, both carry the same x//3 dependence; 
however LLS miss the overall linear dependence on h, accounted for by the function f(2(3h 2 ) in Eq. (8). 
The short distance behavior of Eq. (19) is precisely Eq. (10). 



IV. DEMONSTRATION WITH NUMERICAL SOLUTION 



In this section we demonstrate the validity of Eq. (8) as the continuum limit of the lattice magnetization 

Tr(e^ H / T af) 

(a(j,h B ,T))= Tr(e-^V (20) 

where H is the Hamiltonian of the lattice model, Eq. (1). The quantum boundary Ising chain can be solved by 
applying a Jordan- Wigner transformation, which yields a quadratic fermionic Hamiltonian. The magnetization is 
nonlocal in terms of these fermions: calculation of a(j, h B , T)) is then equivalent to evaluation of the determinant of a 
matrix whose elements are fermionic correlation functions. 20 We construct these analytically, but ultimately evaluate 
them numerically. Details of this calculation follow in Sec. IV A. Here we pre-empt that discussion and present our 
numerical results, comparing to the refined exact expression, Eq. (8). 
The continuum limit expression Eq. (8) admits the scaling form 

(<r(x,h,T)) =T 1 ' 8 J : [x/p,2ph 2 ]. (21) 

For this function to be a continuum limit of the lattice magnetization, there should exist nonuniversal constants 
c, c x , Ch such that 

(<r(j, h B ,T)} = cT^ 8 T[c x j/(3, 2c h ph%] (22) 

is satisfied for all j, hs and T, as long as distances are large compared to the lattice constant, j > 1. and the energy 
scales Kb and T are small compared to the lattice cutoff scale, h Bl T -C L 

The constant c x is related to the velocity v of bulk excitations via c x = v . This follows from the requirement that 
the exponential decay at long distances j 3> f3 is given by 21 (a(x = j)) = (cr(zi = ix)a(z2 = —ix)} — > e^^ 2 ^ 7T ^ 2 ^^ v ^\ 
with v = 1/16 being the scaling dimension of the chiral a field. In our model we obtain c x = 1/2 exactly, c is an 
overall factor relating the lattice magnetization to the field theory one, and Ch relates the (squared) boundary field, 
Kb, in the lattice model to h appearing in the continuum action, Eq. (14). We determine c and Ch by demanding that 
the ratio 

(ajj, h B ,T)) 

cTV^{c xJ /(3,2c h (3h 2 B } [ ' 

is equal to unity for all h B - The best fit from our numerical data was obtained for Ch — 0.161 and c ~ 0.729. 

As shown in Fig. 3, we obtain essentially perfect agreement between numerical calculations and field theoretical 
predictions for the magnetization as a full function of distance, over a wide range of 2h 2 /T = 2chh 2 B /T. 

Fig. 3 should be seen as confirmation that the x dependence of the magnetization is described by the LLS result, 
Eq. (6). However, the full dependence on h, x and T — capturing the evolution from the T <C h 2 result of CZ, 
Eq. (2), to the perturbative T ^> h 2 result, Eq. (19) — is only recovered on inclusion of the factor f(2(3h 2 ) appearing 
in Eq. (8). 
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Figure 3: Comparison of analytical result using Eq. (8) (full lines) and numerical results (dashed lines) for fixed temperature 
T = 0.5 10" 3 and varying boundary field h B = 10" 3+n/4 /%/2, n = 0, 1, 10 increasing from bottom to top [explicitly 
h B = 0.000707, 0.00125, 0.00223, 0.00397, 0.00707, 0.0125, 0.0223, 0.0397, 0.0707, 0.125, 0.223]. 



A. Modified Jordan- Wigner transformation and construction of the magnetization determinant 

We now describe the calculation of the magnetization using a fermionic representation of the transverse field 
quantum Ising chain, Eq. (1). We start from a finite lattice with L sites, 

L-2 L-l 

H l = ~ E «i - E °t - . ( 24 ) 

•i=0 i=0 

with boundary field hs at site j = 0; and with free boundary conditions at site j — L — 1. Ultimately we will take 
the L — > oo limit to avoid finite size effects. 

Consider first the usual Jordan- Wigner representation of the Pauli matrices Tj (j = 0, L — 1), 



fi-i \ 

ilA,ilB,l 1BJ, (25) 



n 

^=0 



in terms of self-Hermitian (Majorana) lattice fermions r yA(B) t ji satisfying {^ajtIaj 1 } = 25^/, {jB,j,lB.j'} — 2<Jjy, 
{7Aj", 7s.j'} = 0. Here, tJ can be obtained from irj = TjTj. Employing this representation for Eq. (24), one obtains 
a linear term involving a single fermionic operator representing the boundary spin operator, Tq = —75,0- This proves 
to be inconvenient in the following, and so we use a modified fermionic representation of the spins to eliminate this 
linear term from the Hamiltonian. Specifically, we introduce an extra boundary Majorana fermion 7 (with = 1), 
which anticommutes with all other fermions -fA.j and Jbj- It can be checked that, if [Tj,T^,] = 2ie abc 5jj'Tj , then 
{crj, (jj, a?} = {a? ,ij<Tj ,iy<Tj} also satisfy [tr",(T^/] = 2ie abc Sjj'<jj. Thus we work with the modified Jordan- Wigner 
representation 

(Tj = ijB,j7A,j, 

\ ijA,ejB,e J ib.]- (26) 
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This is formally equivalent to embedding the spins in a larger Hilbert space. The model Eq. (24) now becomes a tight 
binding model of Majorana fermions, containing only quadratic terms: 

L-2 L-l 

H l = E *7aj7b ,3+1 + E *X4j7l\j + h B ijjB,o- (27) 

3=0 j=0 
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The model can be straightforwardly diagonalized by introducing the fermionic modes 

L-l 

2 



1 

An = o X] C 9 ™^ + Dta,.? + iQntihB.j) , (n = 1, 2, L) 



3=0 



with 



9n(j) = 



L 



■ ( nn , ■ i 



(28) 



(29) 



satisfying the completeness relation £)n=i 9n{j)9n{j') = Z)n=i fl'nO + Dfi'nO'+s) = • This g ives {^n, ^n'} = <S n ,n', 
and 



7A,j 



n=l 
L 

J2g n (j)(A n -Al). 



(30) 



71=1 



The Hamiltonian thus becomes, 



n=l 



n=l 



with E n = 4 cos 2L+i ' which consists of a band of fermionic levels coupled to a Majorana impurity. 
Using Eq. (26) the magnetization is given by 

(of) = -i 3+1 {l TB,o7A,o7B,i7A,i-7B,j-i7-A,j-i7B,j)- 



(31) 



(32) 



One proceeds using Wicks theorem, 20 applicable for the quadratic Hamiltonian Eq. (31). Due to the bipartite structure 
in Eq. (27) it follows that ('jAj^fAj') — {llA.j') = {lB.jlB,j') = 0. All nonzero contractions, including relative signs, 
are then captured by the determinant 



/ i{llBfi) i{llB,i) 



i(lAfl"lB,j) 



\ i(lA,j-llBfl) i{lA,j-llB,l) ■ ■ ■ i(7Aj-l7Bj) / 



(33) 



The calculation of the fermionic correlators in Eq. (33) can be done by exact Green function resummation, treating 
the problem as a noninteracting impurity model. 22 In Eq. (31) we have a quasi-continuum of modes labeled by n 
coupled to a localized impurity state 7. The Green functions for n-modes and for the localized state are defined as 



GnAr) = -a(lt^j)(AV A n ,)), 
G 7 (r) = -(T 7 (t) 7 ), 



(34) 



where 0(t) = e HhT Oe Hlt ) and T is Wick's time-ordering operator. We now construct a perturbative expansion of 

the Green functions in hs, with G(iuj m ) = dTe lu>ml 'G(t) in terms of the Matsubara frequencies u m = nT(l + 2m). 
The zeroth-order Green functions are given by 





n 



G (°)^ rn ) = — . 



(35) 



The full impurity Green function can then be written as G 7 (zu; m ) = [(G7 (iuj rn )) 1 — T,(iuj m )] 1 . Writing the 
boundary term in the Hamiltonian as Hl\i 1b = hs Yln=i 9n{fyl ( A n ) ( 1 J j the exact self energy follows as 



E(«o; m ) = hi J2 9l(0) ( 1 -1 ) G^(iu> m ) ( \) = -h% £ <£(0) 



(36) 



s 



We can now calculate the fermionic correlators entering the determinant Eq. (33). With G(t) = T J^uj e lUmT G{iuj m ), 
the correlators involving the impurity fermion are given by 

*(llB, x ) = »(T7(t = 0+) 7B ,x) = 4 T^e-^ 0+ (jjB, x )u m 

UJm 

= T J2e-^ 0+ Y,9n(x)(l(A n - Al)) Um , (37) 
where we used Eq. (30) in the last equality. Proceeding with first order perturbation theory in /i^, we have 



n=l 



-Th B G 7 (zc m )e— 0+ J2 9n{x)g n (Q) . (38) 



Using the exact expression for G 7 , this becomes 



v^L g„(x)g„(0) 
Z^n=l 

4 i ll B Z^n=l (u m ) 2 +El 



t<ryB,*> = -r/i B ^ e rrTTFi — (A 3 ■ (39) 



1 _l f.2 (g„(0)) 2 - 

Similarly, the {'JaJb) correlators are given by 

*(7a,x7b,x') = -T^e- iw - 0+ £ fln (a: + |)<v(*') ( 1 1 ) 6»n'(*«m) ( _\ ) ■ (40) 

Using standard impurity Green function methods, the full G n n'{ioj m ) Green function is seen to contain two terms, 

G nn ,{iu m ) = 6 nn ,GW(iu m ) + h%&M(iu m ) ( ~* ) . 9ri (0)G 7 ^ m ). 9n ,(0) ( 1 -1 ) G^,^). (41) 
Explicitly, the desired correlator can be expressed as 



i(lA,xlB,x>) = Y^9n{x+ \)g n {x')t&nh ) +2/i|T^]e 

1=1 ^ ' u„ 4 + "'B Z^n" = l (cjJ 1 ) 2 +B ;i 

(42) 

All these expressions are exact for the model Eq. (24) containing two boundaries. We are interested in the effect of the 
boundary j = 0, but not on the boundary at j = L — 1. Thus we proceed by taking the limit L — > oo, which reproduces 
the desired semi-infinite chain. Replacing discrete summations over n by integrals. L +i/ 2 En=i ^ Jo*^ 2 ^ with 
6> = 2^-, Eqs. (39) and (42) become 

iirrsj) = -^E^ t^# oV = 0,1,2,...) 

<<7^7bj0 = <r(? + j' + 1) - - f) + ]T " + , = 0,1,2,...) (43) 



where 



. .. 4 r/ 2 sin(0) sin[(2j + 1)0] 



7T 



n 



cj 2 + (4 cos 6») 2 ' 



/ ■ \ 4 r /2 M 4sin(6»)sin[(2j + 2)6»]cos(6') 



7r J uj 2 + (4cos#) 2 

o-Ci) = -~ f^dO co S [(l + 2i)0]tanh('l^^) 
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The magnetization due to a field Kb applied at the single boundary of a semi- infinite chain at finite temperatures is 
thus given exactly by Eqs. (33) and (43). In practice we evaluate the integrals in Eq. (44) numerically, yielding the 
results presented in Fig. 3. It would be interesting to rederive the field theoretical results by analytic evaluation of 
the determinant Eq. (33) in the continuum limit following the methods of Ref. 5. 



V. CONCLUSIONS 



The crossover physics evinced by the boundary Ising model has been shown to play a key role in a surprisingly diverse 
range of physical problems. 1,8-13 Analysis of the exact universal crossover from free to fixed boundary conditions in 
such problems at finite temperatures thus requires the corresponding scaling functions of the boundary Ising model 
to be known exactly. In this paper we obtained the full scaling function for the magnetization of the boundary Ising 
chain at finite temperature. Among the potential applications of our results, one example is calculation of finite- 
temperature conductance crossovers in two-channel or two-impurity Kondo quantum dot systems. 23 The crossover 
from non-Fermi liquid to Fermi liquid physics in such systems is characterized by the same renormalization group 
flow as occurs in the boundary Ising chain. 10 We plan to extend our earlier work 12 at T — in this area to finite 
temperatures, employing the results of this paper. 23 We note in this regard that without the function f(2/3h 2 ), the 
conductance near the non- Fermi liquid fixed point is unphysical; 23 but using the main result of this paper, Eq. (8), 
exact results 10 at both non-Fermi liquid and Fermi liquid fixed points are recovered precisely. 
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